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Polycomb group (PcG) proteins regulate major developmen- 
tal processes in Arabidopsis. EMBRYONIC FLOWER 2 (EMF2), 
the VEFS domain-containing PcG gene, regulates diverse 
genetic pathways and is required for vegetative development 
and plant survival. Despite widespread EMF2-\\ke sequences 
in plants, little is known about their function other than in 
Arabidopsis and rice. To study the role of EMF2 in broccoli 
(Brassica oleracea var. Italica cv. Elegance) development, we 
identified two broccoli EMF2 (BoEMF2) genes with sequence 
homology to and a similar gene expression pattern to that in 
Arabidopsis (AtEMF2). Reducing their expression in broccoli 
resulted in aberrant phenotypes and gene expression 
patterns. BoEMF2 regulates genes involved in diverse devel- 
opmental and stress programs similar to AtEMF2 in 
Arabidopsis. However, BoEMF2 differs from AtEMF2 in the 
regulation of flower organ identity, cell proliferation and 
elongation, and death-related genes, which may explain the 
distinct phenotypes. The expression of BoEMF2.1 in the 
Arabidopsis em/2 mutant (Rescued em/2) partially rescued 
the mutant phenotype and restored the gene expression pat- 
tern to that of the wild type. Many EMF2-mediated molecu- 
lar and developmental functions are conserved in broccoli 
and Arabidopsis. Furthermore, the restored gene expression 
pattern in Rescued em/2 provides insights into the molecular 
basis of PcG-mediated growth and development. 

Keywords: Brassica oleracea var. italica • Broccoli • 
EMBRYONIC FLOWER 2 • Flowering • Polycomb Group Protein. 

Abbreviations: AG, AGAMOUS; AP3, APETALA3; asBoEMFH, 
broccoli harboring 35S::antisenseBoEMF2.1; CDS, coding sequence; 
CLF, CURLY LEAF; COR75A, COLD REGULATED ISA; EMF2, 
EMBRYONIC FLOWER 2; EXPA, EXPASIN A; FIE, FERTILIZATION 
INDEPENDENT ENDOSPERM; FIS2, FERTILIZATION INDEPEN- 
DENT SEED 2; FLC, FLOWERING LOCUS C, FT, FLOWERING 
LOCUS T; LTP, LIPID TRANSFER PROTEIN; MEA, MEDEA; MS, Mur- 
ahige and Skoog; MSI1, MULTICOPY SUPPRESSOR OF IRA1; PcG, 
polycomb group; PI, PISTILLATA; PRC, polycomb repressive 



complex; qRT-PCR, quantitative real-time reverse transcription- 
PCR; RACE, rapid amplification of cDNA ends; Rescued em/2, em/2 
mutant harboring 3SS::BoEMF2. 7; RFLP, restriction fragment length 
polymorphism; SEP3, SEPATALLA3; Su[z]12, Suppressor of Zeste 
12; SVP, SHORT VEGETATIVE PHASE; SWN, SWINGER; TFL2, 
TERMINAL FLOWER 2; transWT, 35S::BoEMF2.7 containing WT; 
VRN2, VERNALIZATIONS WT, wild type. 

The nucleotide sequences reported in this paper has been sub- 
mitted to GenBank under the accession numbers JQ412719 
(BoEMFZl cDNA), JQ412720 (BoEMF2.2 cDNA), JQ412722 
(BoEMF2.1 genomic DNA) and JQ412723 (BoEMF2.2 genomic 
DNA). The Arabidopsis and broccoli microarray data were sub- 
mitted to Gene Expression Omnibus (GEO) database under the 
accession numbers GSE36943 and GSE36963, respectively. 



Introduction 

Polycomb group (PcG) proteins are epigenetic repressors in 
eukaryotes known to maintain the silent states of their target 
genes (Morey and Helin 2010, Margueron and Reinberg 2011). 
PcG proteins form multiprotein complexes, such as polycomb 
repressive complex 1 and 2 (PRC1 and PRC2), which modify 
histone moieties and remodel chromatin structures to repress 
transcriptional activities (Schuettengruber et al. 2007). PRC2 
has histone methyl transferase activity that trimethylates 
histone 3 lysine 27 (H3K27) on target genes, whereas PRC1 
catalyzes histone 2A (H2A) monoubiquitination (Muller and 
Verrijzer 2009). Together, these PcGs confer stable gene silen- 
cing (Morey and Helin 2010, Margueron and Reinberg 2011). 

PRC2, first described in Drosophila, is composed of four core 
proteins (Muller et al. 2002). In plants, Arabidopsis proteins 
homologous to components of Drosophila PRC2 are required 
for three developmental functions (Calonje and Sung 2006, 
Hennig and Derkacheva 2009, Holec and Berger 2012). 
Arabidopsis FERTILIZATION INDEPENDENT ENDOSPERM 
(FIE) and MULTICOPY SUPPRESSOR OF IRA1 (MSI1) are 
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homologs to Drosophila Extra sex comb (ESC) and p55, respect- 
ively. CURLY LEAF (CLF), MEDEA (MEA) and SWINGER (SWN) 
are homologous of Enhancer of Zeste (E[z]), which catalyzes 
the trimethylation of H3K27 on target genes. EMBRYONIC 
FLOWER 2 (EMF2), VERNALIZATION2 (VRN2) and FERTILIZA- 
TION INDEPENDENT SEED2 (FIS2) are homologous to Suppres- 
sor of Zeste 12 (Su[z]12) (Calonje and Sung 2006, Hennig and 
Derkacheva 2009, Holec and Berger 2012). The three Arabidop- 
sis homologous proteins VRN2, EMF2 and FIS2 share several 
amino acid sequence domains with their Drosophila Su[z]12 
homologs, including the VEFS domain (Yoshida et al. 2001). 
These three proteins also play a major role in determining 
the biological function of the three PRC2 complexes by regulat- 
ing specific target genes. EMF2, together with MSI1, FIE and CLF 
or SWN (EMF2/CLF or SWN/FIE/MSI1), maintains vegetative 
growth by repressing the expression of the flower MADS-box 
genes AGAMOUS (AG), APETALA 3 (AP3) and PISTILLATA (PI) 
(Moon et al. 2003, Calonje et al. 2008). FIS2/MEA/FIE/MSI1 
regulates Arabidopsis seed development by repressing 
PHERES1 during gametophyte and endosperm development 
(Kohler et al. 2003). Likewise, VRN2/CLF or SWN/FIE/MSI1 is 
involved in vernalization-mediated flowering by regulating 
FLOWERING LOCUS C (FLG) (Sung and Amasino 2004, Wood 
et al. 2006, De Lucia et al. 2008). 

Similarly, Arabidopsis proteins homologous to Drosophila 
PRC1 components such as the RING-finger homologs 
AtRING1A/B and AtBMI1A/B are required for maintaining cell 
identity in Arabidopsis (Sanchez-Pulido et al. 2008, Xu and Shen 
2008, Bratzel et al. 2010). These proteins interact with LIKE 
HETEROCHROMATIN PROTEIN 1 (LHP1)/TERMINAL 
FLOWER2 (TFL2) and EMF1, and catalyze histone 2A (H2A) 
monoubiquitination (Bratzel et al. 2010, Chen et al. 2010). 

£A/IF2-like sequences exist in many plants (Chen et al. 2009, 
Luo et al. 2009); however, their function has been studied only 
in the model plants Arabidopsis (Luo et al. 1999, Gendall et al. 
2001, Yoshida et al. 2001) and rice (Luo et al. 2009). Mutations 
in Arabidopsis and rice EMF2 cause early flowering and affect 
flower organ development (Yoshida et al. 2001, Luo et al. 
2009). However, the em/2 mutant is pleiotropic (Yoshida 
et al. 2001) and the mutation affects the expression of genes 
in various genetic programs, including leaf, stem and flower 
development; hormone synthesis and response; and stress re- 
sponses (Moon et al. 2003, Kim et al. 2010). 

To investigate the function of EMF2 in other plant species, we 
studied the close relative of Arabidopsis, broccoli (Brassica oler- 
acea var. italica cv. Elegance) from Brassicaceae. Broccoli is a 
paleopolyploid (Qiu et al. 2009) that diverged from a 
common ancestor with Arabidopsis about 14.5-20.4 million 
years ago (Yang et al. 1999) but differs in gross morphological 
features from Arabidopsis. To study the function of broccoli 
EMF2 (BoEMF2), we cloned two AtEMF2 homologs, BoEMF2.1 
and BoEMF2.2. Both BoE/\AF2 knockdown in broccoli and expres- 
sion of BoEMF2 in Arabidopsis emf2 mutants showed that 
BoEMF2s are similar to the Arabidopsis homolog in being 
required for broccoli growth and differentiation. Unlike 



AtEMF2, BoEMF2s do not seem to promote cell elongation 
and are needed to prevent wilting and premature death. We 
discuss the role of PcG gene modulation in morphological dif- 
ferentiation during evolution. 



Results 




Two BoEMF2 genes identified and 
characterized in broccoli 

Two broccoli (B. oleraceae var. italica) EMF2 (BoEMF2) genes, 
BoEMF2.1 and BoEMF2.2, were cloned by combining 
chromosome walking and rapid amplification of cDNA ends 
(RACE) with primers based on Arabidopsis sequences 
(Supplementary Fig. S1, Supplementary Table SI). 
BoEMF2.1 and BoEMF2.2 are both composed of 22 exons 
(Fig. 1A). Their coding regions are 1,896 and 1,890 nucleotides, 
and show 90.6 and 86% similarity, respectively, to Arabidopsis 
EMF2 (AtEMF2). Southern blot analysis with a probe specific to 
BoEMF2.1 exon 10-14 (Fig. 1A) that showed 81% similarity to 
that of BoEMF2.2 confirmed the cloning of two copies of 
BoEMF2s in the broccoli genome (Fig. 1B). Xba\, Sad, Nco\ 
and EcoRV were chosen to digest broccoli genomic DNA. 
There are two Xba\, two Sad, three Nco\ and one EcoRV cutting 
site on the BoEMF2.1 genome whereas there are three Sad, 
three Nco\ and one EcoRV cutting site on the BoEMF2.2 
genome. As shown in Fig. 1B, digestions of broccoli genomic 
DNA with Xba\, Sad, Nco\ and EcoRV all produced one strong 
and one weak band. The strong band represented the DNA 
fragment of BoEMF2.1 and the weak one represented that of 
BoEMF2.2 (Fig. 1B) due to the fact that the probe is 100% 
match to the BoEMF2.1 sequence. 

Protein sequences of BoEMF2.1 and BoEMF2.2 show 82.2% 
identity to each other, and 90.8 and 83.2%, respectively, to 
AtEMF2. In addition to their high protein identity, BoEMF2s 
and AtEMF2 share a common domain organization (Chen 
et al. 2009): an N-terminal cap, an N-terminal domain, E5-10 
(exon 5-exon10), a C2H2 domain, El 5—17 (exon 15-exon 17) 
and a C-terminal acidic-W/M (VEFS) domain (Fig. 1 C). The VEFS 
domains for BoEMF2.1 and BoEMF2.2 share 87.4% identity to 
each other, and 93.7 and 88.2%, respectively, to AtEMF2. The two 
BoEMF2 proteins are quite different from the other two 
Arabidopsis Su[z]12 homologs, AtVRN2 and AtFIS2 (Chen 
et al. 2009). AtVRN2 and AtFIS2 show only 30.2 and 19.1% pro- 
tein identity to AtEMF2, respectively. Of the protein domains, 
both AtVRN2 and AtFIS2 lost the N-terminal cap region, and the 
E5-10 and E15-17 domains present in AtEMF2 and the 
BoEMF2s. In addition AtFIS2 protein acquired a specific aspartic 
acid-rich (D-rich) domain at the N-terminus and a serine-rich 
(S-rich) domain flanking the C2H2 and VEFS domains (Fig. 1C). 

Real-time quantitative reverse transcription-PCR (qRT-PCR) 
of the expression of the BoEMF2 genes with Brassica Actin 
(BoActin) used as a control revealed ubiquitous expression in 
all organs studied (Fig. 1D). Both genes were most highly ex- 
pressed in flower buds relative to roots, leaves, stems and 
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Fig. 1 Characterization of BoEMF2.1 and BoEMF2.2. (A) Gene structure of the two BoEMF2 genes. The box represents the exon and the 
connecting gray line represents the intron of the two genes. The 5 r - and 3 r -non-coding regions are shown as boxes with diagonal lines. The 
restriction sites of Xba\, Sad, Nco\ and EcoRV, and the region for probe synthesis for Southern blot analysis are indicated. The scale bar represents 
1 kbp of nucleotides. (B) Southern blot analysis. Genomic DNA was digested with use of Xba\, Sad, Nco\ and EcoRV, and detected with a 
BoEMF2. 1 -specific digoxigenin labeling probe on a PVDF membrane. Arrows indicate the enzyme-digested DNA fragments. The size of the marker 
in base pairs is indicated at the left. (C) Protein identity and domain organization. The N-terminal domain (N-ter) is composed of two parts 
including an N-terminal cap (cap), E5-10, C2H2 zinc finger, E15-17 and VEFS domain conserved in At- and BoEMF2s. The D-rich and S-rich 
domain are specific to AtFIS2. Protein identity is indicated at the right. The scale bar represents the length of the protein in amino acids. 
(D) Tissue-specific expression pattern. They-axis indicates the relative mRNA level normalized to BoActin. Data are the mean with SD from three 
independent experiments. 
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Fig. 2 Phenotypes and expression of Bo EMF2 genes in asBoEMF2.1 transgenic plants. (A) Six-week-old vector-only control and three types of 
asBoEMF2.1 transgenic broccoli are shown in the upper panel: type I asBoEMF2.1~41, type II asBoEMF2.1~29 and type III asBoEMF2.1-8. The lower 
panel shows the mRNA level of BoActin and BoEMF2 genes in control and transgenic plants. PCR cycles for BoEMF2.1, BoEMF2.2 and BoActin were 
30, 30 and 24, respectively. The Arabic numbers indicate the expression of genes relative to that of BoActin, used as a control. The arrows 
indicated the curled leaves. (B) Adult plants at the inflorescence stage of control (1), type I asBoEMF2.1 (2) and type II asBoEMF2.1 (3). Elongated 
inflorescence (4), flower bud (6), open flower (7 and 8) of control (left) and type I asBoEMF2.1 (right). An enlarged view of the red box in 2B-4, 
showing a fused pistil-like structure (FP) on the apex of a type I asBoEMF2.1 inflorescence (5). Inflorescence (9) and abnormal flower (10) of type II 
asBo£/V1F2.7. H, broccoli head. P, pistil. Bars represent 1 cm. 



siliques. The expression pattern for both is similar to that of 
AtEMF2. 

Thus, BoEMF2 genes are similar to their Arabidopsis 
homolog in terms of sequence, domain organization and 
gene expression pattern. We knocked down the expression of 
the BoEMF2 genes in broccoli to study their biological function 
and introduced BoEMF2.1 into Arabidopsis em/2 mutants to 
determine whether BoEMF2.1 can rescue the AtEMF2 function. 

Reducing BoEMF2 activity affected vegetative and 
flower development in broccoli 

To investigate the role of BoEMF2 in broccoli growth and 
development, we used reverse genetics to knock down the 



expression of the BoEMF2 genes by introducing the full-length 
coding sequence (CDS) of antisense BoEMF2.1 under the 
control of a constitutive Cauliflower mosaic virus promoter 
35S (35S::antisenseBoEMF2.1) into broccoli. Compared with 
plants transformed with the empty vector (vector-only 
control), T n and T 2 generations of transgenic broccoli harboring 
35S::antisenseBoEMF2.1 (asBoEMF2.1) showed aberrant pheno- 
types during development (Fig. 2; Supplementary Fig. S2). For 
the T n plants, 6 weeks after germination, leaves tended to curl 
upwards. The extent of curl was inversely related to the level of 
endogenous BoE/\AF2 genes (Fig. 2; Supplementary Fig. S2). 
We could classify three types of asBoEMF2. 7 transgenic plants 
by their morphological features. 
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Type I plants showed a few mildly curly leaves and expression 
of BoEMF2.1 reduced as compared with the wild type (WT) level 
(asBoEMF2.1~41 in Fig. 2A; asBoEMF2.P1, -7, -70, -76, -23, -25, -39 
and -47 in Supplementary Fig. S2), whereas the RNA level of 
BoEMF2.2 in these plants was not as affected as in type I plants. 
Type II asBoEMF2.1 plants were dwarf, had tightly curled leaves 
and showed almost completely knocked-out expression of 
BoEMF2.1 and a clear reduction in the asBoEMF2.2 level 
(asBoEMF2.1-29 in Fig. 2A; asBoEMF2.1-14, -29 and -32 in 
Supplementary Fig. S2). Type III asBoEMF2.1 plants had a few 
tightly curled leaves that wilted as new leaves emerged and died 
prematurely within a few days in the soil. Interestingly, unlike the 
emf2 mutation that greatly reduced internode and petiole elong- 
ation in Arabidopsis, type III plants appeared to have a normal 
petiole and internode length (asBoEMF2.1-8 in Fig. 2A). These 
plants had non-detectable BoEMF2.1 levels and almost com- 
pletely knocked-out expression of BoEN\F2.2 (asBoEMF2.1-8 in 
Fig. 2A; asBoEMF2.1-5 and -8 in Supplementary Fig. S2). 

The vector-only control, type I and type II asBoEMF2.1 plants 
could undergo reproductive growth and produce flowers. The 
asBoEMF2.1 broccoli flowered earlier than control plants 
(Table 1): 69.0 ± 5.4 d after germination and producing 
12.8 ±1.4 leaves as compared with 87.7 ± 5.2 d with 16.2 ±1.0 
leaves for controls. Compared with control plants (Fig. 2B-1), 
asBoEMF2.1 type I (Fig. 2B-2) and type II (Fig. 2B-3) plants pro- 
duced thinner stems, fewer leaves, abnormal flowers and, some- 
times, a determinate inflorescence (Fig. 2B-4, -5). These features 
resulted in an abnormal inflorescence profile in asBoEMF2.1 
plants (Fig. 2B-4). 

Type I asBoEMF2.1 plants produced many flowers, thus 
forming a small broccoli head (Fig. 2B-2). The sepals are not 
fully enclosed to protect the inner flower organs in the flower 
bud (Fig. 2B-6). In the open flowers, the petals and sepals were 
irregular in shape and smaller (Fig. 2B-7, -8). Occasionally, 
fused pistil-like structures were observed on the top of the 
inflorescence, which suggests determinate growth 
(Fig. 2B-4, -5). Few seeds could be collected from type I 
plants. Type II asBoEMF2.1 plants produced many to a few 
flower buds (Fig. 2B-9) with imperfect flowers, usually lacking 
petals and stamens (Fig. 2B-10), with barely any seed produced. 

In summary, asBoEMF2.1 showed knocked-down endogen- 
ous BoEMF2.1 and BoEMF2.2 expression, which is consistent 
with the sequence similarity of the two BoEMF2 genes. 



Table 1 Flowering time in asBoEMF2.1 and vector-only control 
transgenic broccoli 



Transgenic plants" 


Days after germination 


Number of leaves 


Vector-only control 


87.7 ± 5.23 


16.2 ±1.03 


asBoEMF2.1 b 


69.0 ± 5.42 


12.8 ±1.42 



0 The number of plants used for analysis was 10 vector-only control and 12 type 

1 (9) and type II (3) asBoEMF2.1 . 

b Type I and type II asBoEMF2.1 broccoli with knockdown expression of en- 
dogenous BoEMF2 genes. Type I plants flowered 67.3 ± 5.29 d after germination 
and produced 13.0 ±1.58 leaves, whereas type II plants flowered after 74±0d 
and produced 12.0±0 leaves. 



Knocking down both BoEMF2 genes caused more severe pheno- 
types than knocking down only BoEMF2.1, so the two BoEMF2 
genes must have redundant functions and have a dosage effect 
on leaf morphologic features, vegetative growth and flower 
development in broccoli. 

Global gene expression pattern of asBoEMF2.1 
transgenic broccoli 

To study the underlying molecular mechanism of the aberrant 
phenotypes caused by the reduced BoEMF2 expression, we 
studied the global gene expression pattern of 11 -day-old 
asBoEMF2.1 and the vector-only control by microarray assay 
(Redman et al. 2004). The 11-day-old asBoEMFH plants with 
one to two true leaves were used to determine BoEMF2 RNA 
levels. Individual seedlings were genotyped to identify those with 
knocked-down expression. Only type I and type II asBoEMF2.1 
broccoli were found but not type III plants because of their 
delayed germination and slow growth. In comparing the gene 
expression pattern for the asBoEMF2.1 transgenic and 
vector-only controls, we identified 2,973 genes up-regulated 
and 3,369 genes down-regulated in asBoEMF2.1 . Similar to the 
Arabidopsis emf2 mutant, asBoEMF2.1 broccoli showed ectopic 
expression of genes in most of the 15 functional categories 
described by Kim et al. (2010) (Table 2; Supplementary 
Table S2). About 16% of genes ectopically expressed in em/2 
were also ectopically expressed in asBoEMF2.1 broccoli, including 
the flower organ identity genes AG and SEPATALLA3 (SEP3), the 
flowering time gene FLOWERING LOCUS T (FT), the seed mat- 
uration LIPID TRANSFER PROTEIN family genes (LTP) and late 
embryo abundant protein genes (LEA), genes encoding transcrip- 
tion factors and genes involved in hormone and stress responses 
(Table 2). The ectopic expression of these flower organ identity 
genes is consistent with the abnormal flower development and 
early flowering of asBoEMF2. 7 broccoli (Fig. 2). 

To confirm the microarray results, we used RT-PCR to exam- 
ine the expression of broccoli SEP3 (BoSEP3), BoAG and BoPI. The 
Affymetrix ATH1 chip experiment showed BoSEP3 and BoAG 
up-regulation, but not BoPI, in asBoEMF2.1 (Supplementary 
Table S2). RNA was extracted from two independent 
11-day-old asBoEMF2.1 transgenic plants (G-i and G 2 ), the 
11-day-old vector-only control (C), and flower (Cf) and silique 
(Csi) from a vector-only control plant. The BoE/\AF2 mRNA levels 
were reduced in the two asBoEMF2.1 broccoli seedlings. The G n 
seedling showed a reduced BoEMF2 mRNA level relative to the 
control: 18% for BoEMF2.1 and 48% for BoEMF2.2, and the G 2 
seedling showed complete knockout of BoEMF2.1 and 35% 
BoEMF2.2 expression as compared with the vector-only control 
(Fig. 3). In the control broccoli, BoSEP3 and BoAG were prefer- 
entially expressed in flower (Cf) and silique (Csi), whereas BoPI 
expression was specific to flowers. In the 1 1 -day-old asBoEMF2.1, 
both BoSEP3 and BoAG were up-regulated but not BoPI, which 
was repressed, as in control seedlings (Fig. 3). The RT-PCR 
results for flower organ identity genes were consistent with 
the microarray results and suggest that BoEMF2.1 and AtEMF2 
regulate similar, but not identical, target genes. 
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Table 2 Categories of genes up-regulated in both em/2 and 
asBoEMF2.1 



Gl G2 C Cf Csi 



Category 



Genes up-regulated by > 2-fold 
in asBoEMF2.1 broccoli 0 



Flower organ identity 
Flowering time 
Seed 

Gibberellin 

Ethylene 

ABA 

Cold 

Stress 

Heat 

Photoreceptor 
Expansin 

Transcription factor 



AG b,c , AP3 b ' c ' d , BELf, SEP3 h ' c 
ATC, FT b ' c ' d 

LEA b,c , LTP3 h ' c , LTP4 b ' c , SAG27 b , 
Protease inhibitor genes b,c 

RGL1 b 

ERF12 b , MBFlC 
HVA22D b ' c , OSTf 
CBFf, DREB1A C 
KINl b ' c , STZ C 



CSN6A C , FUS12 C 
EXPA7 b ' c 
GLS7 b 

COL9 b , MYB96 C , NACl b , NAM b,c 



a See Supplementary Table S2 for additional up-regulated genes. 
b Genes up-regulated in em/2 based on data in the current study. 
c Gene up-regulated in emf2 based on data in Kim et al. (2010). 
d Genes up-regulated only 1.9-fold in asBoEMF2.1 broccoli. 
AF7, NAC DOMAIN CONTAINING PROTEIN 2; AGL57, AGAMOUS-LIKE 57; ATC, 
ARABIDOPSIS THAU AN A CENTRORADIALIS; ANAC087, NAC DOMAIN 
CONTAINING PROTEIN 87; BELT, BELL 7; CBF1, C-REPEAT/DRE BINDING 
FACTOR 7; COL, CONSTANS-LIKE; CSN6A, COP9 SIGNALOSOME COMPLEX 
SUBUNIT 6A, DREB1A, DEHYDRATION RESPONSE ELEMENT B7A; EXP A, 
EXPANSIN A; FUST 2, FUSCA 72; HSFC1, HEAT SHOCK TRANSCRIPTION 
FACTOR CI; HSP, HEAT SHOCK PROTEIN; HVA22D, HVA22 HOMOLOGUE D; 
NAC1, NAC DOMAIN CONTAINING PROTEIN 7; NAM, NO APICAL MERISTEM; 
NIK3, NSP-INTERACTING KINASE 3; MBF1C, MULTIPROTEIN BRIDGING FACTOR 
1C; MYB96, MYB DOMAIN PROTEIN 96; OST1, OPEN STOMATA 7; RGL1, 
RGA-UKE 7; SAG21, SENESCENCE-ASSOCIATED GENE 27; STZ, SALT 
TOLERANCE ZINC FINGER. 



BoEMF2.1 partially rescued the Arabidopsis em/2 
phenotype 

The two BoEMF2 genes are similar to AtEMF2 in sequence, 
protein domain organization and expression pattern. 
Because the expression of BoEMF2.1 was higher than that 
of BoEMF2.2 in broccoli and the former showed higher se- 
quence similarity to AtEMF2 than BoEMF2.2, we used 
BoEMF2.1 to investigate whether it can substitute for the 
AtEMF2 function. We generated a 35S::BoEMF2.1 full-length 
CDS construct and introduced it into Arabidopsis em/2 mu- 
tants to determine whether BoEMF2.1 can restore normal 
seedling and adult development in em/2. Hygromycin- 
resistant transgenic plants harboring the em/2 mutation 
were identified by restriction fragment length polymorph- 
ism (RFLP) genotyping (Supplementary Fig. S3A) as 
described (Calonje et al. 2008). 

Analysis of the phenotypes of the homozygous em/2 
mutants harboring 35S::BoEMF2.1 (Rescued em/2) grown 
under short-day conditions showed that BoEMF2.1 partially 
rescued em/2 phenotypes (Fig. 4, Table 3). 



BoEMF2,l 



8oEMF2.2 



BoSEP3 



BoAG 



BoPI 



BoActin 




1.00 1.00 1.00 1.00 1,00 



Fig. 3 RT-PCR analysis of BoEMF2 genes and flower MADS-box gene 
expression in asBoEMF2.1 and control broccoli. mRNA level of the 
genes in 11-day-old seedlings. and G 2 , two independent 
asBoEMF2.1 transgenic broccoli; C, a vector-only control line; Cf, 
flower of the vector-only control; Csi, silique of the control line. The 
mRNA level was normalized to that of BoActin and is presented as 
Arabic numbers beneath each gel photo. The unsaturated PCR cycles 
for BoEMF2.1, BoEMF2.2, BoSEP3, BoAG, BoPI and BoActin were 33, 33, 
27, 30, 24 and 30, respectively. 



At 7d after germination, em/2 mutants showed short 
hypocotyls, petioleless, and small and oval-shaped cotyledons 
(Chen et al. 1997). BoEMF2.1 partially rescued the hypocotyl 
and petiole elongation in em/2 mutants (Fig. 4A, Table 3). The 
petioleless cotyledons of em/2 were slightly curled upwards, 
whereas the cotyledons of Rescued em/2 were curled down- 
wards along the proximodistal axis, and the WT cotyledons 
were flat and round (Fig. 4A). 

At 15d, WT seedlings produced four petiolated rosette 
leaves and the cotyledons continued to grow (Fig. 4B, C), 
whereas the cotyledons of em/2 stopped growing, and no petio- 
lated rosette leaves were produced (Fig. 4D, E). Rescued em/2 
produced petiolated rosette leaves to varying extents 
(Fig. 4F, G). The hypocotyl was longer for Rescued em/2 than 
em/2 but shorter than for the WT (Fig. 4A-G, Table 3). The 
15-day-old WT and Rescued em/2 were transplanted into the 
soil for further observation, and em/2 seedlings remained on 
agar plates because they cannot survive in soil. The Rescued 
em/2 produced about six petiolated rosette leaves, which were 
smaller than those of the WT (Fig. 4H-K) before flowering at 
about 30 d after germination; em/2 mutants flowered at about 
22 d and produced no rosette leaves, and the WT flowered at 
about 104d and produced >64 rosette leaves (Table 3). 
Therefore, BoEMF2.1 partially rescued vegetative growth in 
em/2, although the rescued plants were smaller and flowered 
early. The Rescued em/2 produced a normal, although smaller 
inflorescence than the WT (Fig. 4J, K). em/2 plants produced 
terminal flowers with a large pistil and dwarfed short stamens 
(Fig. 4L). Rescued em/2 showed indeterminate growth and 
secondary inflorescences developed from the axillary meristems 
(Fig. 4J, K), like WT plants. 
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Fig. 4 Partial rescue of Arabidopsis em/2 by BoEMF2.1. (A) Seven-day-old WT (left), emf2 (middle) and Rescued em/2 (right). (B and C) 
Fifteen-day-old WT. (D) and (E) Fifteen-day-old emf2. (F and G) Fifteen-day-old Rescued emf2. (H) Six-week-old WT at the vegetative stage. 
(I-K) Six-week-old Rescued emf2 at the reproductive stage with varying degrees of rescue in leaf size and vegetative development. The arrowhead 
indicates the emerging inflorescence. The red box in J shows an enlarged view of a petiolated rosette leaf. (L) Six-week-old emf2 produced a pistil 
directly from the embryonic shoot apical meristem (*). (M) A lateral view of the WT flower (left) and Rescued emf2 flower (right) showing the 
smaller petals and shorter stamens. (N) Aerial view of the WT flower (left) and Rescued emf2 flower (right). (O) WT silique (left), with a 
spontaneously developed silique, a short and empty silique in the Rescued emf2 without manual pollination (middle), and a silique of Rescued 
emf2 after manual self-pollination (right). The arrow indicates the cotyledon. The bar represents 2 mm. 



Table 3 Hypocotyl length of 15-day-old Arabidopsis seedlings and 
leaf number at bolting 



Genotype 


Hypocotyl 


Bolting time b 


Rosette leaf 




length 0 


(day, n = 20) 


number at 




(mm, n - 20) 




bolting (n = 20) 


Wild type 


7.0 ± 0.5 


104.4 ±5.4 


64.5 ± 2.2 


em/2 


3.1 ± 0.3 


22.9 ± 2.9 


0 


Rescued em/2 


5.3 ± 0.7 


30.6 ±2.8 


5.9 ±2.5 



a Fifteen-day-old seedlings grown under short-day conditions on agar plates. 
b Arabidopsis transferred from agar at 14 d to soil and grown under short-day 
conditions until bolting. 



Flower development in em/2 was partially rescued by the 
constitutive expression of BoEMF2.1. Rescued em/2 produced 
flower organs of relatively normal size: unlike petal-less and 
extremely dwarfed stamens in em/2, petals were present, al- 
though smaller; the pistil was only slightly taller than stamens 



(Fig. 4M, N). However, partially rescued flowers did not pro- 
duce seeds and the siliques did not elongate spontaneously 
(Fig. 40), probably because of the slightly shorter filaments, 
smaller anthers and fewer pollen produced than the WT. 
Therefore, we performed manual pollination by placing the 
entire anther from the same flower or the same plant onto 
the pistil and generated many seeds in siliques of relatively 
normal size (Fig. 40), but still fewer than the WT. 

Thus, broccoli BoEMF2.1 can substitute for Arabidopsis 
AtEMF2 in performing developmental functions in part at 
every stage of Arabidopsis development. 

BoEMF2.1 partially restored the expression of 
misexpressed genes in the em/2 mutant 

To understand the molecular basis of the partial gene rescue, 
we studied the global gene expression pattern of Rescued em/2 
seedlings. The genotype and expression of BoEMF2.1 were 
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characterized to ensure that the Rescued em/2 plants selected 
based on morphological features were indeed the em/2 mutants 
harboring BoEMF2.1 (Supplementary Fig. S3A, B). RNA was 
extracted from 7~day~old Rescued em/2, em/2, WT, and WT-like 
transgenic seedlings harboring BoEMF2.1 (transWT) for micro- 
array analysis as described in the Materials and Methods. 

We first examined the em/2 expression pattern using the 
22 K Affymetrix ATH1 microarray and found 1,490 and 2,461 
genes up- and down-regulated, respectively, by >2~fold in em/2. 
About 94% of the genes misregulated in the current study were 
detected on the custom GeneChip reported by Kim et al. 
(2010). The two independent data sets showed an overlap of 
393 and 223 genes up- and down-regulated, respectively, in 
em/2 seedlings (Fig. 5A; Supplementary Table S3), which rep- 
resent a robust set of EMF2-regulated genes. As reported by 
Kim et al. (2010), we found misexpression of genes in em/2 in 
the 15 functional categories (Supplementary Tables S3, S4), 
including flower organ identity genes PI, AP3, SEP3 and AG; seed 
maturation genes LTP3 and LTP4; cell growth genes EXPANSIN 
A (EXPA); flowering time genes FT, FLC, SHORT VEGETATIVE 
PHASE (SVP); a COLD REGULATED ISA (COR1SA) gene; and the 
heat-shock protein (HSP) HSP 70 gene. 

We then studied the global gene expression pattern in the 
Rescued em/2. The control used was transWT, which is a trans- 
genic plant harboring BoEMF2.1 in the WT background and 
phenotypically indistinguishable from the WT, and hence a 
better control than the WT, which does not harbor any trans- 
gene. A comparison of the gene expression patterns of Rescued 
em/2 normalized to transWT and em/2 normalized to WT 
revealed 1,122 of the 1,490 genes up-regulated in em/2 (75%) 
no longer up-regulated and 1,899 of the 2,461 down-regulated 



genes in em/2 (77%) no longer down-regulated. Of the robust 
set of EMF2-regulated genes, 257 of the 393 genes up-regulated 
in em/2 (65%) were no longer up-regulated (Fig. 5B) and 158 of 
the 223 down-regulated genes in em/2 (71%) were no longer 
down-regulated (Fig. 5B), which suggests that 35S::BoEMF2.1 
restored the misregulation of many genes in em/2. Further ana- 
lysis showed that the expression of genes misregulated in all 
functional categories was restored, for many, to near-normal 
levels (Table 4; Supplementary Table S4). 

The reduced expression of flower organ identity genes API, 
AP3 and PI, and changed expression of flowering time genes FT 
and SVP are consistent with the normal flower organ develop- 
ment and delayed flowering, respectively, in em/2 harboring 
BoEMF2.1. BoEMF2.1 also restored the expression of the gibber- 
ellic acid synthesis gene GIBBERELLIN 2-OXIDASE 2 (GA20X2), 
stress-related genes C-REPEAT/DRE BINDING FACTOR 1 (CBF1) 
and COR1SA, and transcription factors TEOSINTE BRANCHEDl, 
CYCLOIDEA AND PCF TRANSCRIPTION FACTOR (TCP), NO 
APICAL MERISTEM (NAM) and CUP SHAPED COTYLEDON 
(CUC). Many of these genes are direct targets of EMF1 and 
their expression is repressed by EMF1 and PRC2 (Kim et al. 
2010, Kim et al. 2012). A change in the expression pattern of 
EMF-interacting Arabidopsis genes with BoEMF2.1 suggests the 
ability of broccoli BoEMF2.1 to act as a member of Arabidopsis 
PRC2 in silencing Arabidopsis genes. 

Despite a high percentage of genes with expression restored 
by BoEMF2.1 in em/2, several flower MADS-box genes, AG, PI 
and SEP3, were still up-regulated in Rescued em/2 seedlings 
(Supplementary Table S4). To verify the microarray results, 
we performed qRT-PCR analysis of select genes in transWT, 
WT, em/2 and Rescued em/2. We confirmed the increased 





Genes upregulated in smf2 



Genes down regulated in emf2 




Restoration of upregulated genes (393) Restoration of downregulated genes (223) 

Fig. 5 BoEMF2.1 restoration of genes misexpressed in em/2. (A) Venn diagrams of genes up-regulated (left) and down-regulated (right) by 
>2-fold in em/2 mutants in the current study (yellow) and in Kim et al. (2010) (green). (B) Pie charts of restoration of em/2 misexpressed genes 
by the expression of BoEMF2.1 in em/2 mutants. The number (ratio) of genes with expression restored (blue) and not restored (red). 
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Table 4 Restoration of expression of up- and down-regulated genes in 7-day-old emf2 by BoEMF2.1 



Category 


No. of genes 
investigated 


Genes misex- 
pressed in cm/2 


No. of genes 
misexpressed" 


No. of genes 
with restored 
expression in 
Rescued emf2 h 


Example of genes with expression 
restored to near-normal level 


Flower organ 
identity 


22 


Up 

Down 


13 
0 


11 
0 


PI, SEP1, AP3, SHP7, SEP2, CRC, STK, 
API, SHP2 MBP2, SPL4, CRC 


Flowering time 


20 


Up 


4 


3 


FT, ELF4 






Down 


5 


4 


SVP, AGL37, RAW 


Seed 


108 


Up 


35 


30 


L/r gene family, LtA, ULtuz, aiki, 
SAG21, 2S seed storage gene, GRU2, 
GRU3 






Down 


19 


17 


N/K3, OLE04 


Auxin 


99 


Up 


4 


3 


Auxin-responsive protein genes 






Down 


21 


21 


DFL2, IAA2, IAA30, P/N3, P/N4, P/N7 


Gibberellin 


13 


Up 


3 


2 


GA20X2, RGL2 






Down 


4 


4 


GA-responsive protein genes 


Ethylene 


53 


Up 


7 


7 


ERF12, ethylene-responsive genes 






Down 


9 


9 


EIL1, ERS2, EFE, ERF8 


ABA 


12 


Up 


2 


0 








Down 


2 


2 


ABF1, ABF4 


Cold 


19 


Up 


5 


4 


COR75A 






Down 


3 


3 


CBF7, CBF2, COR15B 


Stress 


46 


Up 


2 


2 


Zinc finger family protein genes 






Down 


6 


6 


Dehydration responsive genes 


Heat 


47 


Ud 


10 


9 


HSP17.6, HSP18.2, HSP70, HSPW1 






Down 


4 


4 


DNAJ, HSFC1 


Photosyntheis 


49 


Up 

r 


0 


0 








Down 


10 


10 


PSII system genes 


Photoreceptor 


28 


Up 


1 


1 


Leucine-rich repeat protein kinase 
gene 






Down 


3 


3 


MYB38/RAX2, ELIP1, PIF3 


Histone 


35 


Up 


1 


1 


histone H7 -3 






Down 


2 


2 


histone H2A 12, histone 2B 


Expansin 


20 


Up 


5 


3 


EXPA7, EXPA18, EXPA20 






Down 


10 


9 


EXPA1, 3, 5, 6, 8, 70, 77, 75, EXPB3 


Transcription 
factor 


542 


Up 


67 


56 


NAC1, AF1, CUC3, COL1, COL9, NAM 






Down 


68 


58 


TCP family genes 



a Genes with >2-fold up- or down-regulated expression in e/n/2. 

b Genes with up- or down-regulated expression in em/2 with expression restored >2-fold by BoEMF2.1. 

ABF, ABSGSIC ACID RESPONSIVE ELEMENT-BINDING FACTOR 7; AF7, NAC DOMAIN CONTAINING PROTEIN 2; AGL23, AGAMOUS-UKE 23; CBF, C-REPEAT/DRE BINDING 
FACTOR; COL, CONSTANS-UKE; COR, COLD-REGULATED; CRC, CRABS CLAW; CRU, CRUCIFERIN; CUC3, CUP SHAPED COTYLEDON3; DFL2, DWARF IN LIGHT 2; DNAJ, 
DNAJ HEAT SHOCK N-TERMINAL DOMAIN-CONTAINING PROTEIN; EFE, ETHYLENE-FORMING ENZYME; EIL1, ETHYLENE-INSENSITIVE3-LIKE 7; ELF4, EARLY FLOWERING 
4; ELIP1, EARLY LIGHT -INDUC ABLE PROTEIN 7; ERF, ERF DOMAIN PROTEIN; ERS2, ETHYLENE RESPONSE SENSOR 2; EXPA, EXPANSIN A; EXPB, EXPANSIN B; GA20X2, 
GIBBERELLIN 2-OXIDASE 2; HSFC1, HEAT SHOCK TRANSCRIPTION FACTOR CI; HSP, HEAT SHOCK PROTEIN; IAA, INDOLE-3- ACETIC ACID INDUCIBLE; MBP2, 
MYROSINASE-BINDING PROTEIN 2; NAC1, NAC DOMAIN CONTAINING PROTEIN 7; NAM, NO APICAL MERISTEM; NIK3, NSP-INTERACTING KINASE 3; MYB38/RAX2, 
MYB DOMAIN PROTEIN 38/REGULATOR OF AXILLARY MERISTEMS 2; OLEO, OLEOSIN; PIF, PHYTOCHROME INTERACTING FACTOR 3; PIN, PIN-FORMED; RAVI, 
RELATED TO ABI3/VP1 7; RGL2, RGA-LIKE 2; SAG21, SENESCENCE-ASSOCIATED GENE 27; SHP, SHATTERPROOF; SPL4, SQUAMOSA PROMOTER BINDING 
PROTEIN-LIKE 4; STK, SEEDSTICK; TCP, TEOSINTE BRANCHED!, CYCLOIDEA AND PCF TRANSCRIPTION FACTOR. 
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expression of flower MADS-box genes AP3 and PI and the seed 
genes LTP3 and LTP4 in em/2 but not AG, and SEP3 restored in 
part with broccoli EMF2.1 (Fig. 6). Similarly, the down-regulated 
cell growth genes EXPA3 and EXPA8 and the flowering time 
gene SVP were restored in part in em/2 with BoEMF2.1. 

In summary, the restoration of misregulated gene expression 
is associated with phenotypic rescue, which suggests that genes 
with restored expression in the Rescued em/2 are at least in part 
responsible for the em/ phenotypes. 



Discussion 

Two BoEMF2 genes show functional conservation 
and redundancy in broccoli 

EA/1F2-like sequences common in plants (Chen et al. 2009, Luo 
et al. 2009) have functions distinct from their derived 
sequences, VRN2 and F/S2 (Chen et al. 2009). The two 
BoEMF2 genes were identified by their high sequence identity 
and their similar gene expression patterns to AtEMF2 (Fig. 1). 
The similar phenotypes between em/2 mutants and asBoEMF2.1 
broccoli (Fig. 2, Table 1) and the partial rescue of em/2 mutants 
by BoEMF2.1 further confirmed their identity. The enhanced 
asBoEMF2.1 phenotypes by knocking down both BoEMF2 
genes implied their functional redundancy. Reducing EMF2 
activities in Arabidopsis and broccoli caused misexpression of 
genes in all 15 functional categories including flower MADS-box, 
hormone signaling and EXPA genes, indicating the conserved 
molecular function of BoEMF2 with that of AtEMF2 in 
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LTP4 EXAP3EXAP8 SVP 



Fig. 6 qRT-PCR analysis of gene expression in Rescued em/2. 
Verification of gene expression based on microarray data by qRT- 
PCR in 7-day-old Arabidopsis seedlings: 3SS::BoEMF2.1 containing 
transgenic WT (transWT), WT, em/2 and Rescued em/2. The y-axis 
represents the relative mRNA level of genes normalized to UBIQUITIN 
10 in log value. Genes analyzed are shown on the x-axis. Data are the 
mean with SD from three independent experiments. 



mediating epigenetic gene silencing (Kim et al. 2010, Kim 
et al. 2012). Luo et al. (2009) reported a T-DNA insertional 
mutant of rice EMF2 that flowered earlier than wild-type rice 
and produced abnormal flowers. Thus, EMF2 is functionally 
conserved in regulating flowering time and organ development 
in Arabidopsis (Yoshida et al. 2001), rice (Luo et al. 2009) and 
broccoli. 

Reducing BoEMF2 activity affected vegetative and 
reproductive development in broccoli 

The relatively weak phenotypes displayed by type I and II 
asBoEMF2.1 broccoli resemble the weak phenotypes of 
Arabidopsis impaired in PcG components, such as elf mutants 
(Schubert et al. 2006), tfl2 mutants (Larsson et al. 1998, Exner 
et al. 2009) and Arabidopsis expressing antisense EMF1 under 
the control of the LEAFY promoter in leaf primordia of 
LFY::asEMF1 transgenic plants (Sanchez et al. 2009). Indeed, 
these plants can produce vegetative leaves but are smaller in 
size, and flower early. In addition, leaves are curly around the 
edges in elf and LFY::asEMFl plants. The curly leaf observed in 
asBoEMF2.1 may be attributed to the ectopic expression of the 
flower MADS-box genes such as AG, because the curly leaf 
phenotype is restored in ag elf double mutants of Arabidopsis 
(Goodrich et al. 1997). The early flowering and abnormal flower 
organ phenotypes of type I and II asBoEMF2.1 broccoli are con- 
sistent with the ectopic expression of flower MADS-box genes 
such as broccoli AG and SEP3 (Fig. 3). SEP3, a class E flower organ 
identity gene, enhances the expression of AP3, PI and AG to 
specify flower organ identity (Immink et al. 2009, Kaufmann 
et al. 2009) and is involved in flowering time control by its 
regulation of flowering time genes such as SUPRESSOR OF 
OVEREXPRESSION OF CONSTANS 7 (SOG7), AGAMOUS LIKE 24 
(AGL24) and SVP in the inflorescence (Kaufmann et al. 2009). 

The microarray data of the type I and type II asBoEMF2.1 
broccoli might display a mild gene expression change due to a 
slight reduction of the activities of BoEMF2s. The strong pheno- 
type displayed by type III transgenic broccoli is, in principle, 
comparable with that of the Arabidopsis em/2 mutant but differs 
from em/2 in that the broccoli seedlings produced 3-5 leaves, 
showed petiole and internode elongation (asBoEMF2.1 -8 in 
Fig. 2A) and died quickly without producing any flower 
organs. Limited availability of type III asBoEMF2.1 plants pre- 
cluded microarray analysis to elucidate the molecular basis of 
these phenotypic differences. However, microarray analysis of 
the type I and type II plants revealed similar as well as divergent 
functions or gene programs in the two species. For example, 13 
flower organ identity genes were misexpressed in the Arabidopsis 
em/2 mutant while eight were misexpressed in asBoEMF2.1 
(Table 4, Supplementary Table S2). Ectopic expression of 
fewer flower organ identity genes may have allowed an abbre- 
viated vegetative development in type III asBoEMF2.1, which 
could not occur in Arabidopsis em/2. Unlike em/2, asBoEMF2.1 
can produce an elongated internode and petiole. This may be 
related to the fact that fewer auxin category and EXPA family 
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genes, which are involved in hypocotyl elongation (Gray et al. 
1998), and cell wall loosening and cell expansion (Sampedro and 
Cosgrove 2005), respectively, are down-regulated in asBoEMF2.1 
(Supplementary Tables S2, S4). In addition, the death-related 
genes, VASCULAR ASSOCIATED DEATH 1 and RADICAL- 
INDUCED CELL DEATH 7, were up-regulated >2-fold in type II 
asBoEMF2.1 broccoli, but remained at a normal level in the 
Arabidopsis emf2 mutant (Supplementary Table S5). This 
may be the reason for the early wilting and death of type III 
asBoEMF2.1 broccoli. 

BoEMF2.1 restored WT seedling and adult 
development in Arabidopsis em/2 

The restoration of nearly 70% genes misexpressed in em/2 
mutants by BoEMF2.1 (Fig. 5) explained the partial rescue of 
hypocotyl, petiole and cotyledon phenotypes of em/2 to that of 
the WT (Fig. 4). These included the restoration of the genes in 
auxin and other hormone categories, as well as most of the 
EXP A genes (Table 4, Supplementary Table S4), which are 
important for cell elongation and expansion in broccoli (Gray 
et al. 1998, Sampedro and Cosgrove 2005). 

BoEMF2.1 also rescued the expression pattern of flowering 
time genes, FT and SVP (Supplementary Tables S3, S4), and 
flower MADS-box genes, thus contributing to the delayed flower- 
ing in the Rescued em/2. Of the 13 ectopically expressed flower 
organ identity genes in em/2, five were no longer ectopically 
expressed and, except for AG and SEP3, the expression of the 
remaining genes was greatly reduced in the Rescued em/2 
(Table 2, Supplementary Table S4). However, Rescued em/2 
still flowers earlier than does the WT, possibly because of the 
ectopic expression of AG and SEP3 (Fig. 6), as SEP3 both activates 
flower organ dentity genes and promotes flowering (Kaufman n 
et al. 2009). The partial phenotypic and molecular rescue of 
em/2 by asBoEMF2.1 may be attributed in part to the use of a 
35S promoter instead of the EMF2 to drive BoEMF2.1 expression. 

Protein-protein interaction between EMF2 and CLF, which 
possesses methyl transferase activity, through the conserved 3' 
VEFS domain in EMF2 is part of the PRC2 multiprotein complex 
that mediates gene silencing in Arabidopsis (Chanvivattana 
et al. 2004). Hence, the fact that BoEMF2.1 restored EMF2's 
function suggests that BoEMF2.1 can interact with the 
Arabidopsis CLF and participate in the Arabidopsis PRC2 
action. This is consistent with its 93.7% sequence identity in 
the VEFS domain to that of AtEMF2. 

However, unlike BoAG and BoSEP3 in broccoli, BoEMF2.1 
cannot repress the expression of AG and SEP3 in Arabidopsis, 
suggesting that these two genes require a more stringent con- 
dition for PRC2-mediated repression than the genes with 
restored expression. As a result, BoEMF2.1 cannot fulfill its 
role in the PRC2 action. The other possibility is that additional 
proteins, such as EMF1, might have evolved in Arabidopsis to 
cooperate with PRC2 in repressing the expression of AG and 
SEP3, but not with the PRC2 containing BoEMF2.1. Future in- 
vestigations at the protein structure-function level will address 



the question of why BoEMF2.1 can restore some but not all 
misexpressed genes in the Rescued em/2, thereby elucidating 
the diverse EMF2-mediated repression mechanisms. 



Materials and Methods 




Plant material and growth conditions 

Seeds of broccoli (B. oleracea var. italica cv. Elegance) were 
purchased from Know-You Seed Company (Kaohsiung, 
Taiwan), germinated on wet filter paper, transferred to soil in 
7 inch pots and grown under a 14/10 h photoperiod at 
75-110 umolrrT 2 s _1 illumination at 20°C. 

Arabidopsis thaliana cv. Columbia (Col) was used as the WT. 
The em/2 mutant constituted the em/2-7 allele in the 
A. thaliana Col background (Chen et al. 1997). Arabidopsis 
seeds were sown directly in soil or surface-sterilized, placed 
on a 0.8% agar plate containing half-strength Murashige and 
Skoog (MS) salt and 1.5% sucrose with or without 20 (igml -1 
hygromycin, cold stratified at 4°C for 2-3 d, then transferred 
to a growth chamber with an 8/16 h photoperiod at 
75-150 jimol m -2 s _1 illumination at 21°C for germination 
and development. Fifteen-day-old seedlings grown on an agar 
plate were transferred into soil and grown under the short-day 
condition for observation of bolting, flowering and silique 
development. 

DNA and RNA extraction 

For Southern blot analysis, genomic DNA was extracted from 
the second leaf from the shoot apex of 1 -month-old broccoli by 
use of the Plant Genomic DNA Extraction Kit (Viogene). For 
gene cloning experiments, genomic DNA was extracted from 
the broccoli inflorescence shoot. 

RNA was extracted from most Arabidopsis and broccoli tissue 
by the TRIzol Reagent method (Invitrogen). To extract RNA from 
siliques, siliques were ground into powder in liquid nitrogen with 
use of a mortar and a pestle, and RNA was extracted with the use 
of extraction buffer (100 mM Tris-HCI, pH 8, 1.4 M NaCI, 20 mM 
EDTA, pH 8 and 2% cetyltrimethylammonium bromide) as 
described (Asif et al. 2000). DNA was eliminated by use of the 
TURBO DNA-/ree kit (Ambion). The RNA remaining was reverse 
transcribed by use of the Transcriptor First Strand cDNA 
Synthesis Kit (Roche) to obtain cDNA products. 

Cloning of broccoli EMF2 (BoEMF2) genes 

Genomic DNA and cDNA prepared from the inflorescence of 
broccoli was used for broccoli EMF2 gene cloning. We used a 
combination of chromosome walking and RACE with the BD 
GenomeWalker Universal Kit (Clontech) and FirstChloice 
RLM-RACE kit (Ambion), respectively. 

Primers designed to clone BoEMF2 genes are given in 
Supplementary Table S1. 

Cloning began with RT-PCR with the primer pairs for EP1 
and EP2 based on the Arabidopsis EMF2 (AtEMF2) sequence to 
obtain a broccoli cDNA fragment of about 500 nucleotides 
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(Supplementary Fig. S1). The following EP primers were de- 
signed with broccoli sequences obtained during cloning. We 
used EP3 and EP4 primers designed with the broccoli 500 bp 
cDNA fragment to obtain the 3' sequence of the cDNA by 
3'-RACE and EP5 and EP6 to obtain the 5' sequence of genomic 
DNA by chromosome walking. With the broccoli genomic DNA 
sequence (the red part in Supplementary Fig. S1 obtained with 
EP5 and EP6), we designed EP7 and EP8 primers to generate 
another cDNA fragment to form a contig with the 3' cDNA 
sequence obtained with primers EP3 and EP4. Further 5'-RACE 
was performed with EP9 and EP10 primers to extend the broc- 
coli EMF2 5'-cDNA sequence. Finally, EP11 containing a start 
codon and EP12 containing a stop codon were designed to 
obtain the full-length coding region of BoEMF2.1. However, 
sequencing results for this full-length BoEMF2.1 coding region 
showed sequence contamination, so we searched for another 
EA/1F2-like sequence (BoEMF2.2) (see below). 

To clone the BoEMF2.2 gene, primers EP13 and EP14 specific 
for a newly identified partial sequence of the £/V1F2-like sequence 
were designed to obtain a cDNA fragment of about 210 nucleo- 
tides (Supplementary Fig. S1). With this short sequence, EP15 
paired with EP16 and EP17 paired with EP18 were designed for 3' 
end and 5' end chromosome walking to obtain the 3' end and 5' 
end genomic sequences, respectively. Further 3'-RACE was per- 
formed with EP13 and EP16 primers to obtain the 3'-cDNA se- 
quence for BoEMF2.2. The following primers were used to 
produce the partially overlapping 5'-cDNA sequences: EP17 
and EP19, based on the BoEN\F2.2 genomic sequence; EP11, 
shared with both BoEMF2 genes; and EP20, specific to the 
BoEN\F2.2 genomic sequence (Supplementary Fig. S1). The 
final 5'-RACE and chromosome walking involved the EP21 and 
EP22 primers to obtain the 5'-untranslated and S'-genomic se- 
quences, respectively. A contig of these cDNA and genomic DNA 
fragment was then reconstructed (Supplementary Fig. S1). 

Southern blot analysis 

In total, 10 ug of total DNA was digested with restriction 
enzymes, separated on a 1.2% agarose gel, blotted onto a nitro- 
cellulose membrane, hybridized with a gene-specific probe and 
imaged by use of X-ray film. The restriction enzymes were se- 
lected with the criterion of no cutting site resides in the gene 
sequences. The gene-specific probes were labeled with digox- 
igenin by PCR (Neuhausurl and Neuhaus 1993). 

Sequence analysis 

All sequences were retrieved from the NCBI (http://www.ncbi 
.nlm.nih.gov/) or TAIR (http://www.arabidopsis.org/) database. 
Sequence analysis involved use of Vector NTI (Invitrogen). 

Plasmid construction, transformation, and 
screening for transgenic broccoli and rescued 
em/2 Arabidopsis seedlings 

Plasmid construction. The Cauliflower mosaic virus promoter 
35S (35S) sequence of pBI121 was amplified by PCR with the 



primers 35S-BaF (TTTGGATCCAGATTAGCCTTTTCAATT) 
and 35S-SaR (AAAGTCGACCGTGTTCTCTCCAAATGA) with 
BamH\ and Sal\ cutting sites, respectively. The amplified 35S 
promoter fragment was then introduced into pCAMBIA1380 
vector by use of the BamYW and Sal\ sites, with selection markers 
of ampicillin for microbes and hygromycin for plants to create a 
pCAMBlA1380'35S'promoter2 vector. 

For BoEMF2 knockdown experiments in broccoli, RNA ex- 
tracted from broccoli inflorescence was used to synthesize 
cDNA. Using the cDNA as template, ABoEMF2-SpF (AAA ACT 
AGTATGCCAGGCATTCCTCTC) and A BoEMF2-SaR (AAAGTC 
GACTCAATTTGGGAGCTGTTC) with Spe\ and Sail cutting 
sites, respectively, were used to clone the full-length CDS of 
antisense BoEMF2.1. The obtained antisense BoEMF2.1 se- 
quence was subcloned into pCAMBIA1380'35S-promoter2 
vector by use of the Spel and Sal\ cutting sites to generate a 
35S promoter-triggered antisense BoEMF2.1 expression vector. 

For the Arabidopsis em/2 rescue experiment, primer pairs for 
BoEM F2-SaF2 ( AA AGTCGACGGGCG ATTATGCCAGGCATTCC 
TCTC) and BoEMF2-SpR (AAAACTAGTTCAATTTGGGAGCTG 
TTC) with Sal\ and Spel cutting sites, respectively, were used to 
clone the sense BoEMF2.1 gene from broccoli inflorescence 
cDNA. The cloned BoEMF2.1 sequence was then subcloned 
into pCAMBIA1380'35S'promoter2 vector by use of the Sal\ 
and Spel cutting sites to create a 35S promoter-triggered 
BoEMF2.1 expression vector. 

Transformation of antisense BoEMF2. 7 into broccoli. Cotyledons 
and hypocotyls of 3-day-old broccoli seedlings were cut into 
pieces and infected with the Agrobacterium tumefaciens strain 
GV3101 harboring 35S::antisenseBoEMF2.1 as described (Chen 
et al. 2001 ). Transgenic broccoli was screened by ability to grow 
on 20|igmr 1 hygromycin and plasmid specific-PCR. 

emf2 rescue by BoEMF2.1. Because homozygous em/2 is lethal, 
we introduced 35S::BoEA/IF2.7 into heterozygous em/2 plants, 
which were identified by allele-specific RFLP genotyping 
(Calonje et al. 2008) and transformed with A. tumefaciens 
strain GV3101 harboring 35S::BoEMF2.7 by the floral dipping 
method (Clough and Bent 1998). T n seeds from transformed 
heterozygous em/2 plants (T 0 ) were surface-sterilized and 
screened on an agar plate containing 20 jig ml -1 hygromycin. 
Hygromycin-resistant seedlings were transplanted into soil to 
obtain T 2 seeds and plants. T 2 plants homozygous for 
35S::BoEA/1F2.7 and harboring the Arabidopsis em/2 mutation 
were identified by producing 100% hygromycin-resistant 
progeny and allele-specific RFLP genotyping. 

RT-PCR analysis 

The second leaves of 4-week-old antisense BoEMF2.1 broccoli 
(asBoE/\AF2.1) and vector-only control were harvested, and the 
expression of the two BoE/\AF2 genes was determined to con- 
firm the knockdown efficiency of the endogenous BoEMF2 
genes (Fig. 3). In total, 2 jig of total RNA was used to synthesize 
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the first-strand cDNA, and PCR analysis involved primers 
specific to BoEMF2.1 and BoEMF2.2 (Supplementary Table 
S1). The expression of broccoli actin (BoActin; accession No. 
AF111812) was a control. 

To study gene expression in asBoEMF2.1 broccoli, 11-day-old 
seedlings of the asBoEMF2.1 T n generation (G-,), the T 2 gener- 
ation (G 2 ), vector-only control seedlings and flower and siliques 
of control plants (Fig. 4) were harvested for RNA extraction and 
cDNA preparation. In total, 2 ug of total RNA was used for 
first-strand cDNA synthesis. The cDNA was diluted five times 
in contraction before PCR analysis with gene-specific primers 
(Supplementary Table S1). The PCR products were separated 
on a 1.5% agarose gel and stained with ethidium bromide. 
Images were obtained by use of a Bio-Rad Doc EQ system, 
and the intensity of DNA bands was quantified with use of 
Quantify One (Bio-Rad). 



qRT-CR analysis 

For qRT-PCR analysis of BoEMF2 genes in broccoli, roots and 
second expanded leaf and stem were harvested from 4-week-old 
broccoli grown in soil; flower buds and siliques were collected 
after flowering and silique formation. For Arabidopsis, RNA was 
extracted from 7-day-old seedlings of the 35S::BoEMF2. 7 -contain- 
ing WT (transWT), WT, em/2 and Rescued em/2. 

Primers for genes were designed for qRT-PCR (Supple- 
mentary Table S1) with the ABI Prism 7500 sequence detection 
system (Applied Biosystems) at 58°C for 2 min, 95°C for 10 min, 
and 40 cycles of 95°C for 15 s and 60°C for 1 min. The relative 
gene expression was determined by the comparative CT method, 
with the expression of UBIQUITINW (UBQW, At4g05320) used as 
an internal control for Arabidopsis and BoActin for broccoli. The 
mean value of 2 _ACT (AC T = C T , gene of interest - C T , control 
gene) was determined from three independent experiments. 

Microarray experiments and analysis 

For broccoli, the seedlings of asBoEMF2.1 transgenic plants and 
vector-containing control were grown on half-strength MS agar 
containing 1.5% sucrose and harvested at 11 d after germin- 
ation. RNA was extracted by the TRIzol reagent method 
(Invitrogen), and RT-PCR was performed to determine the 
endogenous BoEMF2 level. The asBoEMF2.1 line with 
knocked-down expression of BoEN\F2s and the vector-only con- 
trol (Fig. 3) were used for microarray analysis. 

For Arabidopsis, 7-day-old 35S::Bo£/VlF2.7-containing WT 
(transWT) and Rescued em/2 seedlings grown on hygromycin- 
containing agar plates were harvested and sorted by cotyledon 
phenotype: transWT and Rescued em/2 with round flat and 
curled outwards cotyledons, respectively (Fig. 4). The seedlings 
of WT and em/2 grown on MS-only plates were collected and 
sorted by round flat and small oval shaped cotyledons, respect- 
ively (Fig. 4). After confirming the genotype by allele-specific 
RFLP (Calonje et al. 2008) and BoEMF2.1 expression 
(Supplementary Fig. S2), microarray analysis was performed. 



For microarray experiments, 10 jig of total RNA was used as 
a template for cDNA synthesis, cRNA was biotin labeled by 
in vitro transcription then fragmentation by use of the 
GeneGhip Expression Analysis Technical Manual rev5 
(Affymetrix). In total, 10 ug of labeled cRNA was hybridized 
to the 22 K ATH1 GeneChip (Redman et al. 2004) at 45°C 
for 16.5 h. Washing and staining involved use of Fluidic 
Station-450, and genechips were scanned by use of the 
Affymetrix GeneChip Scanner 7G. The mean gene expression 
indices for all chips were scaled to an arbitrary target of 500. 

For broccoli gene expression analysis, only genes with hy- 
bridization signals >60 in at least one of the two samples, the 
vector-only and the asBoEMF2.1 broccoli, were analyzed. For 
Arabidopsis, only genes with hybridization signals >50 in at 
least of one of four plant samples, em/2, WT, Rescued em/2 
and transWT, were analyzed. The hybridization signal ratio of 
fold change of gene expression in broccoli was analyzed by the 
ratio of asBoEMF2.1 to vector-only control, with a ratio >2 or < 
0.5 used to identify up- or down-regulated genes, respectively. 
The fold change of gene expression in Arabidopsis was based on 
the ratio of em/2 to the WT over that of Rescued em/2 to 
transWT, with a ratio >2 or <0.5 used to identify up- or 
down-regulated genes, respectively. 



Supplementary data 

Supplementary data are available at PCP online. 
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